Hysteresis in a system driven by either generalized force or displacement variables 
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We report on experiments aimed at comparing the hysteretic response of a Cu-Zn-Al single crystal 
undergoing a martensitic transition under strain-driven and stress-driven conditions. Strain-driven 
experiments were performed using a conventional tensile machine while a special device was designed 
to perform stress-driven experiments. Significant differences in the hysteresis loops were found. The 
strain-driven curves show re-entrant behaviour (yield point) which is not observed in the stress- 
driven case. The dissipated energy in the stress-driven curves is larger than in the strain-driven 
ones. Results from recently proposed models qualitatively agree with experiments. 

PACS numbers: 81.30.Kf,75.60.Ej 
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I. INTRODUCTION 

Hysteresis is ubiquitous in many areas of physics. It is 
a signature of non-equilibrium effects and has attracted 
the attention of researchers for a long time and is still one 
of the most challenging subjects of research^. Typically, 
the phenomenon of hysteresis results in closed loops when 
plotting the external driving force (or field) vs. the gener- 
alized displacement (or vice versa). These two variables 
(external force and generalized displacement) are conju- 
gate variables in the first principle of thermodynamics 
and depending upon the specific system, they can have 
different tensorial character. Nevertheless, in all cases 
their product renders the work performed on the studied 
system and the area enclosed within a closed loop gives 
the energy loss. 

Among hysteretic systems, particular attention has 
been devoted to athermal cases^. In these cases, hys- 
teresis is not due to competition between the fast driv- 
ing rate and slow thermal relaxation, but rather it orig- 
inates from the existence of very high energy barriers 
that can only be overcome when the system reaches lo- 
cal marginal stability limits^i^. The continuous change 
in the driving force results in discontinuous changes in 
the conjugate displacement, giving rise to avalanche dy- 
namics. Such dynamics has been observed in a wide va- 
riety of materials such as ferromagnetic^, ferroelectric^ 
and martensitic- materials; field-driven vortex motion in 
type-II superconductors^, and pressure-driven conden- 
sation of '^He on mesophorous solids^, among others. 
Since most of these systems show common hysteretic be- 
haviour they are catalogued under the name of fcrroic 
materialsi^. They display first order phase transitions 
with an order parameter (exhibiting discontinuous be- 
haviour), which corresponds to the generalized displace- 
ment. A common feature in these systems is the ex- 
istence of long range effects arising from compatibility 
conditions. 

From an experimental point of view, it is usually easier 
to control the force while the conjugate displacement is 
measured. For instance, in magnetic systems, the mag- 



netic field is readily controlled and magnetic flux is mea- 
sured. Although more difficult, it has also been possi- 
ble to control the magnetic flux^^ - which is suitable for 
studying materials displaying a rapid increase of mag- 
netization. Moreover, the control of the amount of gas 
(generalized displacement) instead of the chemical poten- 
tial has been shown to be the most convenient procedure 
in a number of gas adsorption experiments^^. 

For a macroscopic system in equilibrium, thermody- 
namic trajectories do not depend on which is the control 
variable (generalized force or displacement) because the 
two cases are related by a Legendre transformationi^. 
However, there is no reason to expect that out of equi- 
librium hysteresis loops obtained when controlling the 
force will be similar to those obtained when the control 
parameter is the generalized displacement. Such a situa- 
tion has been very recently theoretically studied^^ii^ by 
making use of the Random Field Ising Model (RFIM)^^, 
which is a prototype model for athermal hysteresis. Re- 
sults predict significant differences in the hysteresis loops 
obtained in the force-driven case (external field H) to 
those obtained by driving the generalized displacement 
(magnetization M). It is the aim of the present paper to 
address this issue from an experimental point of view. 

Martensitic transitions (MT) offer a unique scenario to 
undertake such a task. A MT is a displacive first-order 
transition which involves a change in symmetryii. In 
these ferroic systems, the transition can be induced by 
cooling but also by application of a mechanical stress. In 
this latter case, the driving force is the applied load, while 
the conjugate displacement is the elongation (or strain). 
Typical experiments are carried out using commercial 
tensile testing machines in which the control variable is 
the elongationi^. We have developed an experimental de- 
vice which enables fine control of the applied force while 
the strain is monitored. In this paper we present experi- 
ments performed on the same specimen under both force 
and displacement control conditions. Results will enable 
a meaningful comparison of the hysteresis loops obtained 
in the two cases. In addition, comparison with the pre- 
dictions of recent theoretical approache a^'^d^ will also be 
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presented. 



II. EXPERIMENTAL 



A Cu-Zn-Al single crystal was grown by the Bridg- 
man technique. The nominal composition was chosen so 
that the (athermal) MT on cooling without stress takes 
place from a cubic {L2i) to a monoclinic (18i?) structure 
slightly below room temperature (Tm = 234K). The ac- 
tual composition, obtained from electron dispersion anal- 
ysis is Cu68.13Zn15.74Ali6.13. A sample was mechanically 
machined from the ingot with cylindrical heads. The 
body of the sample has flat faces 35 mm long, 1.4 mm 
thick and 3.95 mm wide. The axis of the sample is close 
to the [001] crystallographic direction of the cubic phase. 
The sample was mechanically polished and then annealed 
for 20 min at 1073K, cooled in air down to room temper- 
ature and aged for 2 h in boiling water. This heat treat- 
ment ensures that the sample is in the ordered state, free 
from internal stresses and that the vacancy concentration 
is minimum at room temperature. 

Special grips, which adapt to the heads of the speci- 
men, can be used in both experimental devices. In the 
first device [Fig. [^a)], a INSTRON 4302 tensile machine, 
(from now on strain- driven), the control parameter is the 
elongation while the load is continuously monitored. A 
second device (from now on stress- driven) was especially 
designed which enables control of the load applied to the 
sample while the elongation is continuously monitored 
[see Fig. [Hb)]. The device was adapted from that pre- 
viously used by Carrillo et ali^. The upper grip is at- 
tached to a load cell hanging from the ceiling. Upper 
and lower ball-and-socket joints ensure good alignment. 
The lower grip holds a container that plays the role of a 
dead load. The load can be increased or decreased at a 
well-controlled rate by supplying or removing water by 
means of a pump. The strain is measured by a strain 
gauge attached to the sample. Since the stress needed 
to induce the MT is very sensitive to temperature, it is 
important to make sure that equivalent experiments are 
carried out at the same temperature. For this reason we 
use a cryofurnace (which enables temperature control to 
an accuracy of 0.1 K.) that can be adapted to both de- 
vices. In addition, in order to be able to compare the 
hysteresis loops obtained with the two devices, the same 
strain gauge was used in all experiments and the load cell 
for each set-up was calibrated using standard weights. 
For selected experiments at room temperature (without 
the cryofurnace), an optical microscope was incorporated 
into both devices, which enables in situ observations of 
microstructural changes. 



III. RESULTS AND DISCUSSION 

Fig. shows typical results for the stress-strain curves 
obtained with the two devices. Upon increasing load. 
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FIG. 1: Schematic representation of the strain-driven (a) and 
stress-driven (b) experimental devices 
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FIG. 2; (Color online) Experimental stress-strain hysteresis 
loop at 303.1 K in a CussZnieAlig single crystal, obtained 
in strain-driven (squares) experiments at 0.005 mms~^ and 
stress-driven (circles) experiments at 0.4 Ns~^. The inset 
shows an enlarged view of the low-strain region. 



the parent cubic phase persists until it becomes unsta- 
ble at a given value of the load, and the MT starts. 
Different hysteresis loops are observed in the two cases. 
Almost the entire strain-driven loop is enclosed within 
the stress-driven loop, which has a much larger area. In 
the stress-driven case, the stress slightly increases dur- 
ing the MT with an average slope that decreases as the 
driving rate is reduced. Extrapolation to zero rate still 
gives a finite (small) slope which can be caused by a 
weak concentration gradient along the axis of the sam- 
ple. Different behaviour is observed in the strain-driven 
experiment. In this case, once the transition starts, the 
stress relaxes to a lower value (yield point effect) and 
then weakly oscillates around an almost constant value 
as the MT progresses. For both driving modes, when 
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FIG. 3: Micrographs showing the microstructure at different 
stages of the MT: (a) 0.45% , (b) 0.48 % and (c) 0.57% strain. 
They illustrate the re-entrant behaviour in the strain-driven 
experiment, (d) 101.6 MPa, corresponds to the stress-driven 
case and illustrates nucleation of thin martensitic plates. 



the transition ends, a new reversible elastic regime in 
the martensitic phase is reached. Upon unloading, the 
behaviour parallels that observed in the loading branch, 
but with hysteresis. In particular for the strain-driven ex- 
periment, before the end of the reverse transition, there 
is an increase in the stress which reproduces the yield 
stress effect. 

Optical observations have revealed significant differ- 
ences in the early stages of the MT as shown in Fig. [3l 
The inset in Fig.[2]presents an enlarged view of the upper 
branch of the hysteresis loops on the low-strain region. 
The yield point in the strain- driven experiment is a con- 
sequence of re-entrant behaviour within the parent phase, 
associated with shrinkage and eventually the disappear- 
ance of previously formed martensitic plates. Such re- 
entrant behaviour is also observed on unloading, which 
causes the yield point of the lower branch of the hystere- 
sis loop. The evolution of the microstructure associated 
with this re-entrant behaviour is illustrated by the micro- 
graphs in Fig. [3](a),(b) and (c). The existence of a yield 
point in conventional strain-driven experiments has also 
been reported for other martensitic alloys^Si^. More- 
over, computer simulations of the mechanical response 
of a strain-driven martensitic system also show a yield 
poin1)2^. 

For the stress-driven experiment there is a deviation 
of the pure elastic behaviour well below the yield stress. 
This effect is due to the nucleation of a set of very thin 
martensitic plates all over the sample, as illustrated by 
the micrograph in Fig. [^d). Interestingly, this is the 
only region where the strain-driven loop is outside the 
stress-driven one. Actually, when driving the load, nucle- 
ation can occur at lower values since the system can more 
freely adapt to shape changes. For both driving mecha- 
nisms, steady state growth following the early stages is 
macroscopically similar, with parallel interfaces growing 
towards both ends of the specimen. 

It is worth mentioning that the behaviour found in 



the present experiments for MT seems to be rather com- 
mon in other systems. For instance, magnetization- 
driven magnetic hysteresis loops also show re-entrant 
behaviour^^ while no re-entrancy is present in magnetic 
field-driven loops. Furthermore, recent theoretical ap- 
proaches and computer simulation studies for athermal 
plastic deformation'^^ provide curves with yield points 
when the control variable is the strain. Interestingly, wig- 
gly trajectories and re-entrance have even been reported 
to occur in displacement driven nanoscale systems such 
as deformation of gold nanowire^ and RNA unfolding^-. 

A relevant quantity in hysteretic systems is dissipated 
energy which is given by the area enclosed within the 
hysteresis loops. A noteworthy feature (see Fig. [2]) is 
that the dissipated energy is much larger for the stress- 
driven case than for the strain-driven one. The dissipated 
energy is expected to depend on the driving rate of the 
control parameter. We have performed a scries of exper- 
iments at different rates and the results of the area of 
the loops are plotted in Fig. |3] as a function of the rate. 
The dissipated energy increases with increasing rate in 
both cases. Extrapolation to zero rate with least-squares 
fits (lines) indicates non-vanishing hysteresis. Similar be- 
haviour has also been reported for magnetic systems^^. It 
is also important to note that even in the zero-rate limit, 
the dissipated energy for the stress-driven experiment is 
always larger than for the strain-driven case. 
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FIG. 4: (Color online) Area of the hysteresis loops as a func- 
tion of the driving rate for the stress-driven case (circles) and 
the strain-driven case (squares). The inset shows the corre- 
sponding results obtained by simulation of the T=0 RFIM 

We now discuss the experimental results presented here 
in relation to two recent approaches (in magnetic termi- 
nology) aimed at analysing _ff-driven and Af -driven tra- 
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jectories in the 3d-RFIM. Both models predict the exis- 
tence of a yield point in the Af -driven loop in agreement 
with experiments. In the iJ-driven case the force cannot 
relax and therefore the system can only reach the closer 
metastable states by a discontinuity in the gencralitzed 
displacement. In contrast, for the Af -driven case the gen- 
eralized displacement is constrained but the force is able 
to relax to lower values in such a way that the system is 
able to follow a trajectory that approaches equilibrium, 
thus resulting in a lower dissipation. The mechanism to 
relax the force is expected to depend on the details of the 
specific dynamic features of each system. The present 
experiments show that in MT the relaxation is achieved 
by retransformation. This is also the mechanism in the 
model in Ref.™. The wiggly trajectories (see Fig. [2]) 
can be interpreted within the framework of the models 
as being due to the specific distribution of disorder in 
the system. One of the models^^ uses T ~ adiabatic 
(infinitely slow) dynamics for the Af-driven case which is 
the analogue of standard adiabatic dynamics introduced 
by Sethna et alfi^ in the 7f-driven case. The definition 
of the measured (output) field assumes that the exter- 
nal force instantaneously equals the highest value of the 
internal forces. Consequently, the Af-driven loops have 
vanishing area. The strain-driven experiments presented 
here indicate that although the area is quite small com- 
pared with the force-driven situation, after extrapolation 
to zero driving rate, a certain amount of dissipated en- 
ergy still remains. 

We have extended the T = RFIM to incorporate the 
effect of a finite driving rate. For the ff-driven case such 
an extension was proposed^l by assuming an intrinsic av- 
erage time for the avalanches Tav to relax. In this case the 
finite-rate ff-driven loops can be constructed from the 
adiabatic ff-driven loops by considering the discontinu- 
ous step-like behaviour of the driving field with a certain 
Aff in such a way that ff = AH /Tav This method al- 
lows the area of the loops as a function of the driving rate 
to be computed. Results obtained for a cubic system with 
size L = 20 and standard deviation of the random-fields 
a = 1.0 are shown in the inset of Fig. |4] (circles). Data 
correspond to averages over more than 1000 realizations 
of the random fields. We have considered a similar defi- 
nition in the simulations of finite-rate Af-driven systems. 
It is assumed that the force is only able to relax at given 
values of the magnetization which are separated by Am 
so that m = Am /Tav An example of the obtained results 
for the same system is also shown in the inset of Fig. 0] 
(squares). Results from numerical simulations satisfacto- 
rily reproduce the major features found experimentally, 
i.e. an increase in dissipation with increasing rate, and 
the fact that the dissipated energy is always larger for 
the force-driven cases. Only at low-rate limit do numeri- 
cal data show the already mentioned tendency towards a 
vanishing area. Such a tendency is due to the assumption 
of the infinitely fast response of the system assumed in 
the models and therefore it is not corroborated by exper- 
imental data. The finite-rate models exactly predict that 



dissipation in the large-rate limit of the Af-driven dynam- 
ics reaches the zero-rate (adiabatic) limit of the ff-driven 
case (see inset of Fig. S]). Experiments should provide a 
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FIG. 5: (Color online) Detail of the rate dependence of the 
strain driven loops at high rates, as indicated. 

test of this prediction. As shown in Fig. [5l experimental 
data indicate a tendency to avoid re-entrance as rate is 
increased (the yield point is less and less pronounced) in 
such a way that the corresponding strain-driven trajec- 
tories approach the zero-rate stress-driven one. This is in 
agreement with numerical simulation predictions. How- 
ever, it must be taken into account that release (or ab- 
sorption) of latent heat (which needs to be driven away) 
associated with the first-order nature of the transition, 
also contributes to dissipation as evidenced by the in- 
creasing slope with increasing rate in the strain-driven 
curves (see Fig. [SJ . Actually, this extra effect cannot be 
reproduced by the model which does not take latent heat 
into account. 



IV. CONCLUSION 

We have measured the stress-strain curves during the 
martensitic transition of a Cu-Zn-Al single crystal. Ex- 
periments have been conducted under well controlled 
force and displacement conditions by means of specifi- 
cally adapted experimental set-ups. This enables a com- 
parison of the hysteresis loops obtained under different 
control variables. It has been shown that metastable tra- 
jectories strongly depend on the driving mechanism. In- 
terestingly, the displacement driven loops exhibit yield 
point and lower dissipation. In situ optical microscopy 
has shown that the mechanism for re-entrance is the re- 
transformation of part of the material. Results essentially 
conform to recent theoretical predictions for athermal 
systems. In spite of the fact that the microscopic mecha- 
nism determining the features of the hysteresis loops can 
be specific to each particular system, it is expected that 
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the main results reported in the present study will be 
common to many ferroic hysteretic systems undergoing 
athermal first-order phase transitions. 
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